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gluino and its maximal number of signal events are about 11, although we cannot realize such model 
in the known supersymmetry breaking scenarios. 
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I. INTRODUCTION 

With the discovery of a 125 GeV Standard Model (SM)-like Higgs boson [Il|2], all the SM particles have already 
been found, and the main goal of the LHC experiment switches to search for physics beyond the SM (BSM). As we 
know, supersymmetry (SUSY) is still the best-motivated framework. It provides an elegant solution to gauge hierarchy 
problem, achieves gauge coupling unification at a high scale, and naturally has the lightest supersymmetric particle 
(LSP) as a dark matter candidate if R parity is preserved. However, the current null results of SUSY searches from 
the LHC together with the measured 125 GeV Higgs mass have placed stringent exclusion limits on the parameter 
space of its simplest realization: the Minimal Supersymmetric Standard Model (MSSM) [3]. Thus, it is worth to 
consider seriously the possible extensions of MSSM, such as the Next-to-Minimal Supersymmetric Standard Model 
(NMSSM) and its extensions. 

Meanwhile, several recently observed excess of signal events may have already given us the new physics clues. 
Notably, the ATLAS Collaboration reported a 3cr excess in the channel of two same-flavour opposite-sign dileptons 
with an invariant mass around Z boson 81 GeV < m£+g- < 101 GeV, jets, and missing transverse energy [!]• The 
goal of this search is to probe two scenarios 

• The gluino g pair productions and the following two-step decays through sleptons to the LSP neutralino x)*, 
which leads to the off-Z dileptons. 

• The general gauge mediation (GGM), where gluinos first via 3-body decay into quark pairs and neutralino Xij 
and the latter then decay into a very light gravitino G plus a Z boson, which results in on-Z dileptons. 

At the 8 TeV LHC with an integrated luminosity of 20.3 fb~^, the observed number of leptonic-Z events for combined 
electron and muon channels is 29 while 10.6±3.2 events are expected in the SM |4], which corresponds to 3cr deviation. 
On the other hand, the CMS Collaboration has also implemented similar search and no significant signal on Z 
excess was observed. The maximal deviation was found 2.6cr in the dilepton mass window 20GeV < < 70GeV. 

Several solutions to explain this excess for both SUSY [6HT5] and non-SUSY models [16] have been proposed. Most 
of SUSY interpretations are based on gluino/squark pair production processes, which can be summarized as follows 

• The MSSM with GGM |7|, the Z-signal can be produced via decay chain g qqXi —> qqZG, where the first 
step decay through an off-shell squark exchange. 

• The MSSM with a light right-handed sbottom bi, a bino-like LSP Xi and nearly degenerated higgsino-like next- 
to-LSP (NLSPs) X 2 3 [Hj- In this scenario, corresponding decay chain is 5 —>■ bb\ —>■ &&X 2 3 bbxfZ. However, 
such rich b signals are tightly constrained by b-jet searches. 

• The MSSM with a tunned spectrum, i.e., light squraks q ~ 500 — 750 GeV, bino-like NLSP X 3 ~ 350 GeV and 

higgsino-like LSPs Xi ,2 ~ 150 — 200 GeV. The corresponding decay chain is g —)■ qxs —t 9 X 1 , 2 -^ El- 

• The MSSM with a split spectrum, the gluinos decay into higgsino-like NLSPs through t-ti-loop, followed by 

higgsinos decay into the bino LSP plus Z boson, i.e., g —>■ 5 X 2 3 9XiZ [TO] . 

• The MSSM with mixed stops El, the Z boson can be produced by heavy stop decay, <2 Zii , where the light 
stop has a mass close to the LSP in order to evade the LHC search constraints. 

• The goldstini model with gauge mediated SUSY breaking scenario in which they also considered gluino 

radiative decay into higgsino-like neutralino g —>■ 9 X 1 , 2 - The difference is there exists a pseudo-goldstino G' in 
the spectrum which leads to Xi ,2 ZG'. 

• Finally, the NMSSM with a singlino-like LSP Xi and a bino-like NLSP X 2 [EHEj) where extra Z bosons can 

be produced through g —>■ 99 X 2 99-^Xi or 9 —>■ 9 x 2 9Xi^- 

In the explanation using light squarks mils], the Majorana gluino mass was chosen by hand to be heavy at low 
energy scale. However, we cannot realize such scenario in the UV completion of the MSSM/NMSSM for the known 
supersymmetry breaking mechanisms since the squark masses will be driven to the order of gluino mass due to the 
renormalization group equation running (For an example, see Ref. [17j.L In this work, we investigate an extension of 
the NMSSM with a Dirac gluino, which provides a UV completion of the above scenario. In particular, the heavy Dirac 
gluino will not induce the SUSY electroweak fine-tuning problem, and then such extension is still natural [I8l - l58] . We 
demonstrate the leptonic-.Z excess can be successfully explained through light squark pair productions. 

Before studying the model in details, we would like to discuss the decay patterns of the first two generation squarks 
in the context of the MSSM and explain our strategy to address the Z excess. With heavy gluino, squarks decay 
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into quarks plus neutralino/chargino: q —> qx^ or g —)■ q'xf- The decay channel q —>■ qxi is always kinematically 
favored and for right-handed squarks it can be dominant since in most case Xi is bino-like. However, the left-handed 
squarks strongly prefer to decay into the wino-like charginos/neutralinos due to large squark-quark-wino couplings. 
Moreover, squarks decay into higgsino-like charginos/neutralinos are only important for the third generation squarks 
in which squark-quark-higgsino couplings are sizeable because of large Yukawa couplings. Due to these properties, 
when the LSP neutralino is bino-like, the NLSP neutralino is wino-like (higgsino-like), squarks can not have large 
branching ratio for cascade decay q —)■ 9X2(X2 3 ) QXi^- While in case of the bino-like NLSP and higgsino-like LSP, 
one may potentially has considerable branching ratio for above cascade decays but the specific fine-tuning of the mass 
spectrum is required. Unfortunately, above situation is also true for the MSSM with a Dirac gluino. So we pay our 
attention to the NMSSM with the Dirac gluino. In this case the LSP (NLSP) neutralino is singlino-like (bino-like) 
and both q —?> qx 2 and X 2 Xi^ '^an have large branching ratios in large parameter space. 

This work is organized as follows. In Section [iTj we present our model and demonstrate its key features. In 
Section III we describe the cutfLow of the ATLAS search for the leptonic-Z-|-jets-|-.^T signal and give the corresponding 
best fit benchmark points. We scan the relevant SUSY parameter space to study the capability of our model in 
interpreting the ATLAS Z-peaked excess by considering the constraints from both the ATLAS and CMS SUSY 
searches. As a comparison, we also discuss the same signal channel in the usual NMSSM within the same parameter 
space. Finally, our conclusion is given in Section [Tvl 


II. THE NMSSM EXTENSION WITH A HEAVY DIRAC GLUINO AND LIGHT SQUARKS 


Dirac gauginos have been proposed decades ago [TBH^ . and recently have been studied extensively in model 
building and phenomenology pTHKH] . The SUSY models with Dirac gluinos have notable advantages comparing 
with their Majorana counterparts. Firstly, they have special renormalization properties, the so-called supersoft 
behavior [HHH sn m] which allows specific mass spectra where squarks are much lighter than glunio. Secondly, 
some subprocesses involved in the first two generations of squarks pair productions vanish (those final states with 
squarks of the same chirality) since the f-channel Dirac gluino exchange cannot mediate the chirality-flip processes. 
This leads to the reduction of squark pair production cross sections in Dirac case compared to the Majorana case, 
which will significantly alleviate their bounds from the current LHC direct SUSY searches, especially in the heavier 
Dirac gluino region {Mg >2 — 3 TeV) [HI [HI SHI SH]- This is known as supersafe behavior and play a crucial role 
when using our model to explain the ATLAS Z-peaked excess. In terms of supersafe, the first and second generations 
of squarks can be as light as 700 GeV, which actually removes the constraints as we want. Therefore, the appropriate 
and convenient framework for weak scale SUSY is the unification of Dirac gluinos and semi-soft SUSY breaking. In 
order to generate the Dirac gluino mass, we need to introduce a chiral superfield with SU{3)c x SU{2)l x U{1)y 
quantum number ( 8 , 1 , 0 ). 

In this paper, we consider the general NMSSM framework whose superpotential is given by [59] 

W = VuQffuU - YdQHdd- YeLHae + (^ + A5)i7„ • Ha + ^fS + (1) 


where Q, u, d, L, e, and S are the NMSSM chiral superfields, A and k are the dimensionless Yukawa couplings, /r 
and are the bilinear mass terms, and is the tadpole term. The corresponding SUSY breaking soft terms are 


= T^QH^u - TdQHj- T^LH^e + ■ H^S + + (gS + h.c. 


m 




\H„, 














+ MiBB + M 2 WW + M 3 GG -k M^BS + M‘gG^ -k h.c., 


( 2 ) 


where Tu,d,e,x,K are the trilinear soft terms, Bg,m'g are the bilinear soft terms, Mi^ 2,3 are respectively the soft 
Majorana gaugino masses for bino, wino and gluino, as well as the Dirac type of gluino (bino) mass {M^). 
Noticed that in Eq. [^ one could has mixed mass gluino eigenstates, i.e., pseudo-Dirac gluino, while the pure Dirac 
or Majorana gluino corresponds to M 3 = 0 or M'g = 0. For simplicity, we consider the pure Dirac gluino and set 
M3 = M^= 0. 

In Fig.^ we show the total cross sections of the first two generation squark pair productions {qq + qq* +q*q*) at the 
8 TeV LHC. For comparison, the NLO-kNLL cross section of Majorana gluino pair production based on NLL-f ast [60] 
calculation is also presented. We find that the cross sections of total squark pair productions are larger than 10 fb 
for Trig < 900 GeV. Without taking into account cut efficiency and branching ratio suppression in each step of squark 
cascade decay chain, we may obtain sufficient Z events in this mass region. The crucial point is then how to keep the 
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Z yield in squark cascade decays and escape the constraints from the ATLAS and CMS direct SUSY searches. To 
achieve this goal, we consider the following strategy 

• The neutralino sector is similar to that of Refs. mE], i.e., the bino-like NLSP X 2 singlino-like LSP x°. 
By choosing M 2 ^ /r 3> Mi, the branching ratios of squark decays into wino and Higgsino are suppressed. As a 
consequence, the branching ratio of squark decay into bino is almost 100%. 

• For simplicity, we assume slepton i and the third generation squarks ti, 2 /^i ,2 are decoupled. Also, to keep the 
nice feature of heavy Dirac gluino model, gluino mass is fixed at rrig > 2.6 TeV. 

• The mass splitting between X 2 and x? is taken to be m^o — m^o ~ 100 GeV to forbid X 2 decay, which 

guarantees Br(x 2 —-^Xi) ~ 100% 0 

In the next section, we present our simulation results for both Z excess signal events and LHC constraints. 



M (GeV) 


FIG. 1: The total NLO+NLL cross section of squark pair productions aysLO+NLL = + QQ* + q*q*) in our model 

at the 8 TeV LHC (red dashed line). For comparison, the NLO+NLL cross section of gluino pair production for 

decoupled squarks is also shown (black solid line). 


III. THE Z-PEAKED EXCESS AND LHC SUSY SEARCH CONSTRAINTS 

For the purpose of studying collider phenomenology, we implement the model in the Mathematica package SARAH[6T]. 
which generates SPheno[5^ output file to calculate the particle mass spectra, mixing matrices as well as the low 
energy constraints. It can also create the corresponding UFO [53] model file which can be used by MG5 _eiMC [M]. In our 
simulation, the parton level events are generated by MG5_aMC, whose output LHE file then feeds into PythiaS [65] to 
implement particle decays, parton showering and hadronization and the package NLL-fast is used to calculate the 
NLO+NLL fc-factor of the squark pair productions. Moreover, we use DelphesS m for fast detector simulation. 

The ATLAS and CMS Collaborations have carried out copious SUSY searches at the LHC. The null signals for most 
of the searches so far have put quite stringent bounds on many SUSY models. As a result, those scenarios, which could 
potentially explain the Z boson excess, may have already been excluded. In order to use the LHC SUSY searches to 
constrain our interesting models, we use a private package in which most of the current LHC SUSY searches have 
been recasted. The package was first used in Ref. and further developed in Refs. iMiizni- The check of validation 
can be found therein as well. 


^ This can also be realized by choosing specific SUSY parameters such that the coupling hxiX 2 suppressed. 

^ In order to employ PythiaS to correctly handle the decays of all SUSY particles in our model, we use the so-called “semi-internal 
processes” which was introduced in Ref. ESj. 
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Specifying to our current work, we find that the LHC searches for the final states with jets and fr and with 0 [7T] 
and without [721I73] isolated leptons are related to the searches for the final states with on-shell Z boson [74]. We 
recast those searches by strictly following the analyses in the corresponding conference notes or published papers. 
The Delphes with default ATLAS setup is used for fast detector simulation. 

In order to measure whether a model is excluded by the current LHC searches or not, we define variables i?* = 

where the N^p and N^p are the numbers of new physics event and upper limit in each signal region Si for a given 
analysis, respectively. The largest = maxi(i?®) among all signal regions of all relevant searches is specified to 

the corresponding model. As a result, R!^^^ > 1 means that there is at least one signal region in which the number 
of signal events is larger than the upper limit of current data, so this model is excluded. In the LHC analyses, the 
ATLAS Collaboration has explicitly given A^p for each signal region, while the CMS Collaboration only presented 
the observed numbers of events and corresponding background events and its uncertainty. We calculate the 95% C.L. 
Npp for the CMS analyses via standard Bayesian procedure as follows 

Nup 

Cin,ts\N,,Nt,ab)P{N,)dN, = 0.95, (3) 

where P{Ns) is uniform prior probability and Af = C{nohs\Ns,Ni,,ab)P{Ns)dNs is a normalization factor. As¬ 
suming Gaussian distribution for background and signal uncertainties, the likelihood function in Eq. [^can be written 
as 



PiRobs I Ag, A^, 


1 

v/2^y2^ 




dbPijiohs] 


(4) 


Practically, the probability function P{nobs'^^J) is taken as Possion distribution ^/nobs^- for Uobs ^ 100 while 

taken as Gaussian distribution for Uobs > 100. Here, the expectation value /i = Ag-|-ds-|-A;,-|- 

and the signal error as = 0.01 rig. We find the derived Ayp from Eq. ® will not change much as long as the signal 
uncertainty is within a few tens of percent. The calculated A^pS for all signal regions in Ref. [5| are shown in Table 
for illustration. 


TABLE I: The 95% C.L. App for all signal regions in the CMS search for events with two leptons, jets, and missing 

transverse momentum |5]. 


Ajeta 

> 2 

> 3 

(GeV) 

100 - 200 

200 - 300 

> 300 

100 - 200 

200 - 300 

> 300 

Total background 

1204 ± 106 

74.5 ± 11.3 

12.8 ±4.3 

478 ± 43 

39.2 ±6.6 

5.3 ±2.3 

Data 

1187 

65 

7 

490 

35 

6 

Aup 

211 

23.2 

8.6 

107 

16.5 

8.9 


With these necessary tools, we are able to apply relevant LHC SUSY search constraints on all benchmark models. 
We now briefly describe the cuts flow of the ATLAS analysis on Z excess. According to the search, the signal events 
should satisfy the following requirements 

1. Events contain at least two same-flavour opposite-sign leptons. 

2. The leading lepton > 25 GeV and sub-leading lepton > 10 GeV. 

3. The invariant mass of these two leptons must fall into Z boson mass window 81 < mi+i- < 101 GeV. 

4. Events further have > 2 jets with pip > 35 GeV and |? 7 | < 2.5. 

5. The missing transverse energy Ifir > 225 GeV. 

6. Ht > 600 GeV with Ht = Pp^ + Pt^ + YhiPr ■ 

7. A(j){ji^ 2 ,$T) > 0.4, where Acj) is the azimuthal angle between two objects in parenthesis. 

In Table [H] we present the most relevant information for two benchmark points. Their cut efficiency, survival 
number of signal events in the ATLAS cut flow as well as the corresponding R^^^ values are shown in Table |Tll| From 
this table, we find that benchl can give the best fit signal events (18.7) and satisfy all of the ATLAS constraints 
we have considered, but not the GMS constraints. While bench2 can escape the constraints form both the ATLAS 
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and CMS searches. As the price, the signal events drop to 14.2, which roughly corresponds to 2cr deviation from the 
SM and is within Itr region of the observed excess. Finally, Fig. shows the Ht and fx distributions respectively 
for two benchmark points after applying all cut selections in the ATLAS cut flow except the Hx and Ipx cuts. Some 
important features can be learnt from here 


From Table III we find that Hx and Ipx cuts are crucial in the cut flow, and hard Hx and fx distributions are 
required in order to keep sufficient survival events. 




As shown from Table [III and Fig. bench2 has relatively compressed spectrum which leads to softer Hx and 
fx distributions compared with benchl. 


As a consequence, bench2 is easier to evade the constraints from the LHC searches especially when the CMS limits 
are taken into account. Meanwhile, fewer events can pass Hx and Ifix cuts in the ATLAS leptonic-Z + jets + ^x 
search. On the other hand, benchl easily satisfies the requirement of Z-signal events but is challenged by the CMS 
constraints. 


TABLE II: The relevant particle masses, branching ratios, and NLO+NLL cross sections for two benchmark points 

of our model. 



M~1,2 

^L,R 

M^o 

too 

Mh 

Br((7L,fl qX 2 ) 

Br(x^ ^ Z ^ i ) 

cnlo+nll (fb) 

benchl 

693 

283 

398 

125 

100% 

95.7% 

79 

bench2 

642 

313 

428 

125 

100% 

95.6% 

133 


TABLE III: The cut efficiency and signal events for our model in the ATLAS cut flow at the 8 TeV LHC with 
luminosity C = 20.3 fb“^. The values are also presented for the ATLAS searches alone as well as for both the 

ATLAS and CMS searches. 



All Events 

Cutl 

Cut2 

Cut3 

Cut4 

Cuts 

Cut6 

Cut7 

^Sig 

Tpmax 

-^ATLAS 

Tpmax 

-^ATLAS+CMS 

benchl 

I(P 

3869 

3839 

3549 

3482 

1599 

1252 

1165 

18.7 

0.56 

1.26 

bench2 

10^^ 

3851 

3822 

3530 

3441 

996 

574 

528 

14.2 

0.52 

0.96 




MET 


EIG. 2: Hx (left) and fx (right) distributions after implementing the ATLAS cut flow, except the cuts on Hx and 
fx- The benchl and bench2 are plotted by black and green solid lines, respectively. 

We further implement a grid scan in the rriq — m^o plane. Contours of signal events and corresponding exclusion 
limits are shown in Fig. Among all searches undertaken in this study, we find that the CMS search for final states 
with jets, ^x and two opposite-sign same-flavor leptons at the Z pole gives the most stringent bound on most of 
the models that could explain the ATLAS excess because of the similar signature they are looking for. The exclusion 
curves in Fig. are shown with (blue curve marked as ATLAS-I-CMS) and without (green curve marked as ATLAS) 
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the CMS constraint, respectively. From the figure, we conclude that after taking into account all the current LHC 
SUSY searches, our model can give at most ~ 15 events for the ATLAS Z-peaked excess, which means the excess can 
be addressed within Icr level. 



m~q [GeV] 


FIG. 3: Contours of signal events for the ATLAS leptonic-Z + jets + Ipr search in the rriq — m^o plane (orange 
dashed curves) for our model. The corresponding exclusion limits for the ATLAS constraints alone (green curve) 
and the ATLAS+CMS constraints (blue curve) are shown as well. 


It is interesting to investigate the conventional low energy NMSSM with heavy Majarona gluino, and compare it 
with our scenario, although such NMSSM cannot be realized in the usual SUSY breaking scenarios. For this purpose, 
we use the same signal channels and mass spectra as in Dirac gluino case. The corresponding results are presented in 
Fig.El Form this figure, we learn the following facts. Firstly, these contours for numbers of signal events have similar 
shape, but the absolute values are distinct. It is mainly due to the larger production rates of squark productions in 
conventional NMSSM. Because the cut efficiencies has no huge discrepancy for these two models, the survival event 
numbers are increased significantly. Secondly, for the same reason, the exclusion limits become stringent and thus 
requiring more compressed spectra. To be specific, for given G [600, 900] GeV, m^o should heavier than 520 GeV 
if one only includes the ATLAS constraints (see green curve in Fig. |^. Within these mass regions, the conventional 
NMSSM can also gives the best explanation for leptonic-Z excess, and the corresponding best fit benchmark point 
is given in Table |IV| as benchS. However, the GMS constraints become so stringent that almost entire parameter 
space in the plane is excluded. Although a few points do survive, the y con tributes about 11 signal events in the most 
optimistic case (About 1.2 (t signal significance, see bench4 in Table [lV| ), and their numbers are too small to draw 
the exclusion curve. Therefore, only the ATLAS constraints are displayed in Fig. The results are similar with 
Ref. |15| . in which the wider squark and neutralino mass ranges were studied. 


IV. CONCLUSION 

In the NMSSM extension with a heavy Dirac gluino and light squarks which can still keep the naturalness condition, 
we interpreted the 3cr leptonic-Z excess recently reported by the ATLAS Collaboration. We concentrated on light first 
two generation squark pair productions, and the corresponding decay chain is g —>■ gy® QXi^- Our basic strategy 
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TABLE IV: The conventional NMSSM with the same caption as Table 



Af-1,2 

^L,R 

M^o 

too 

All Events 

Cutl 

Cut2 

Cut3 

Cut4 

Cut5 

Cut6 

Cut7 

A'sig 

r)max 

-^ATLAS 

omax 

-^ATLAS+CMS 

benchS 

738 

489 

602 

10*' 

14652 

14533 

13462 

12876 

2746 

1314 

1219 

18.5 

0.84 

1.53 

bench4 

665 

511 

624 

n? 

11941 

11843 

10959 

7916 

874 

398 

349 

10.6 

0.48 

0.95 



650 700 750 800 


m~q [GeV] 


FIG. 4: The conventional NMSSM with the same caption as Fig. Noticed that the ATLAS+CMS limits are too 

stringent to be displayed in the figure. 


to produce .Z boson is using a singlino-like LSP Xi and a bino-like NLSP With specific sparticle spectra and mass 
splitting TO^o — m^o ~ 100 GeV, the branching ratios of decays q —)■ qx 2 and x\Z are almost fixed to 100% in 
the whole parameter space. Meanwhile, to satisfy the LHG SUSY search constraints and maintain the nice supersafe 
feature of Dirac gluino, a relatively heavy glunio is chosen, i.e., rrig ~ 2.6 TeV. 

Gonsidering the bounds from all ATLAS direct SUSY searches, the leptonic-Z excess can be fully addressed in large 
parameter space of our model, i.e., G [650,750] GeV and m^o G [200,400] GeV. The benchl with 2=639 
GeV and TOj^o=398 GeV can provide the 19 signal events for the leptonic-Z excess while satisfies all the current 
ATLAS SUSY search constraints (7?“tlas = 0.56). 

The corresponding GMS search for exactly the same final states imposes very stringent bound on the above viable 
parameter space which can explain the leptonic-Z excess. The parameter space, which satisfies both the ATLAS and 
CMS constraints, can provide at most 15 events to the excess, which means the excess can be addressed within la. 
The bench2 with rriq.^ 2=0^2 GeV and mj^o=428 GeV, which contributes 15 signal events to the excess, is on the edge 
of the CMS exclusion curve. 

For comparison, we considered the conventional low energy NMSSM with a heavy Majorana gluino and light squarks 
in the same parameter space. The main difference is the enhanced squark pair production cross sections from the 
same chirality squark pair. In the same mass regions considered in this work, we found the capacity to explain the 
leptonic-Z excess in both models are similar. If one only considered the ATLAS SUSY searches, a quite large viable 
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parameter space can give ~ 19 signal events to the Z excess, e.g., benchS. While the constraints from the CMS 
search only allow the parameter space with less than ~ 11 signal events, e.g., bench4. 
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